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ABSTRACT 

Direct imaging searches have begun to detect planetary and brown dwarf companions and to place 
constraints on the presence of giant planets at large separations from their host star. This work 
helps to motivate such planet searches by predicting a population of young giant planets that could 
be detectable by direct imaging campaigns. Both the classical core accretion and the gravitational 
instability model for planet formation are hard-pressed to form long-period planets in situ. Here, 
we show that dynamical instabilities among planetary systems that originally formed multiple giant 
planets much closer to the host star could produce a population of giant planets at large (ss 10 2 AU 
— 10 5 AU) separations. We estimate the limits within which these planets may survive, quantify the 
efficiency of gravitational scattering into both stable and unstable wide orbits, and demonstrate that 
population analyses must take into account the age of the system. We predict that planet scattering 
creates a population of detectable giant planets on wide orbits that decreases in number on timescales 
of ~ 10 Myr. We demonstrate that several members of such populations should be detectable with 
current technology, quantify the prospects for future instruments, and suggest how they could place 
interesting constraints on planet formation models. 

Subject headings: planetary systems: formation — imaging — techniques: high angular resolution — 
celestial mechanics — methods: n-body simulations — open clusters and associa- 
tions: general 



1. INTRODUCTION 

Recently, several high-contrast imaging campaigns 
have begun to search for young giant planets orbiting 
relatively nearby stars at large orbital separations. In 
addition to compelling detect i ons of likely planeta ry sys- 
tems (e.g. iKalas et all 120081: iMarois et all 120081) . sur- 
veys of young stars have placed significant constraints 
on t he presence of any giant planets at large separations 
("e.g. iMcCarthv fc Zuckermanl l2004t i Kasper etaTl l2007t 
ILafreniere et al.ll2007t iNielsen et atl l2007t) . In the near 
future, even more powerful searches have the potential 
to detect even lower mass planets in wide orbits around 
nearby stars. These surveys could reveal new popula- 
tions of extrasolar planets, since planets with long (> 10 2 
yr) periods are difficult to detect by most planet search 
techniques such as radial velocities, astrometry, transit 
photometry, and even microlensing. 

The frequency and orbital properties of long-period 
planets could provide interesting constraints for planet 
formation models. The core accretion model strug- 
gles even to explain the in situ accretion of Uranus 
and Neptune, given the need to form a rocky 
core while there is still significant gas in the cir- 
cums olar nebula ([Pollack et al.1 Il996t iGoldreich et al.1 
2004). This timescale problem has led to sev- 
eral models in which Uranus and Neptune form 
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giant planets at much larger separations such as ~ 1000 
AU or even ~ 100 AU. Alternatively, the gravitational in- 
stability model may become more plausible at large sep - 
arations ()Rafikovil2005tlStamatellos fc Whitworth| [2008'l. 
Hence, some recent ground-based direct detection cam- 
paigns have framed their searches as a way to search for 
planets that could not be formed by core accretion, and 
thus provide support for some planets having formed via 
gravitational instability. However, recent research has 
also cast doubt on the viability o f the gravitational in- 
stability mode l (iMeiia et alj l2005t ICai et alj 120051: IBossI 
120061 : iRafikovl [2007m . Therefore, one might worry that 
a non-detection by direct imaging planet searches (e.g. 
| Kasper et alj 120071 : ILafreniere et al.l 120071 : INielsen et al.l 
l2007h would fail to provide new constraints on planet 
formation models. In this paper, we show that this is 
not the case. In either the core accretion or gravita- 
tional instability model, the formation of multiple giant 
planets at moderate separations (~ 1 — 10 AU) naturally 
leads to dynamical instabilities that result in one or more 
planets being scattered into very wide (« 10 2 AU — 10 5 
AU) orbits. Thus, direct im aging searches (such as GPI 
The Gemini Planet Imager, Macinto sh et al.ll2006a| ) and 
microlensing planet searches have the potential to place 
significant constraints on the frequency of dynamically 
active planetary systems that form multiple giant plan- 
ets. 

Extrasolar giant planets have a broad range of 
eccentricities. Although several proposed mecha- 
nisms may explain the eccentricities of some planets, 
the high frequency of eccentric orbits among gi- 
ant planets suggests that dynamical instabilities 
resulting in planet-planet scattering play an im- 
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portant role in sculpting many planetary systems. 
Although many variations of the plane t-planet scat- 
terin g model of been explored (e.g . iRasio fc Fordl 
19961: IWeidenschilling fe Marzaril Il996t iLevison et al.l 
19981: iFord etalJjggOlt iMarzari fc Weidenschillind 120021: 
Adams fc Laughlinl 120031 : IVeras fe Armitaeel 1200 ' 



2006: iMoorhead fc Adamsf 120051: IChatteriee et all 1200 



Ford fc Rasiol l2008t Uuric fc Tremaind 12008ft . these 
studies focused on the final properties of the bound 
planets and paid little attention to the orbital evolution 
of planets that will eventually be ejected from the 
system. Recently, a dynamical instability has been 
suggested as a possible formation mechanism for the 
subst ellar companion to GQ Lup (jDebes fc Sigurdssonl 
2006) in a wide (« 100 AU) orbit. However, previous 
dynamical studies have yet to quantify the probability 
for scattered planets to be detected via direct imaging. 
We investigate the orbital properties and lifetime 
of planets scattered into wide orbits and assess the 
prospects for their detection vi a high-contrast imaging . 
Our study complements that of IScharf fc Men ou (2008), 
as we consider additional effects, including galactic tides, 
scattering of terrestrial-mass planets, and observational 
estimates of the contrast for long-period giant planets. 

Given the history of extrasolar planets repeatedly be- 
ing discovered in unexpected locations, it is worth search- 
ing for planets anywhere that they could survive. There- 
fore, we begin by exploring the physical limits that can 
be placed on very long-period planets and find that plan- 
ets could survive with orbital separations of up to ~ 10 5 
AU in ^5] Such estimates can be compared to the sep- 
arations at which planets are observed in order to de- 
termine whether there are additional processes shaping 
this distribution. In <j3l we describe a model where mul- 
tiple giant planets form at traditional orbital separa- 
tions (^3 — 30 AU). As the disk dissipates, the planets 
will perturb each other, eventually resulting in eccen- 
tricity growth, close encounters, and some planets be- 
ing ej ected from the planetary system ([Chatteriee et al.l 
2008). Before the planets are ejected, they will spend 
considerable time at very large orbital separations from 
the host star. Thus, planet scattering should result in a 
significant number of giant planets being scattered onto 
very wide orbits, where they could be detected by re- 
cent or future direct imaging planet searches. By quan- 
tifing the efficiency of scattering onto such orbits, we 
estimate how large a sample needs to be surveyed in or- 
der to detect such planets. We additionally explore the 
timescale on which detectable wide-orbit planets remain; 
we predict that the population of unstable planets "pass- 
ing through" a wide-orbit phase is larger than the popu- 
lation of planets on stable wide orbits. The existence of 
this population does not require that either the standard 
core accretion or gravitational instability mechanism be 
able to form giant planets at large separations, where di- 
rect imaging campaigns are most sensitive. Observations 
of these planets at different ages can impart information 
about the chaotic stage of planet formation (e.g., how 
many giant planets are typically formed before planet 
scattering widdles down the number of surviving planets 
to 2 or 3). Finally, we consider the prospects for a direct 
imaging campaign to detect such planets and their impli- 
cations in 21 System age, contrast, and separation limits 
all affect population analyses from such campaigns. 



2. SURVIVAL AND DESTRUCTION OF LONG 
PERIOD PLANETS 

We begin by asking what processes would limit the 
survival of very long-period giant planets orbiting main- 
sequence stars. Most previous studies investigating the 
effects of passing st ars have focused on planets at separa- 
tions < 40 AU (e.g.lLaughlin fc Adamdl 19981 : IFord et all 
l2000at I Adams et al.l 1200a ) . If giant planets have much 
larger orbits, then the effects of stellar encounters could 
be dramatically increased. In the limit that planets are 
much less massive than stars, they can be treated as test 
particles, allowing us to apply many arguments previ- 
ously applied to the survival of long-period comets. 

2.1. Field Stars 

First, we consider the constraints on long period plan- 
ets around field star s in the ga l actic disk, closely follow- 
ing the treatment of iTrem ainc (1993|). 

2.1.1. Galactic Tide 

Planets in sufficiently wide orbits will be influenced by 
gravitational perturbations from the galactic tides. Like 
comets, a planet will typically become unbound if the 
star-planet separation exceeds the tidal radius, estimated 

by 
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where M* is the mass of the host star , and p is the local 
galactic mass density ()Tremainelll993| ). Once a planet is 
scattered beyond a t , it will likely become unbound on an 
orbital timescale. 

2.1.2. Passing Stars 

Encounters with passing stars and giant molecu- 
lar clouds can eject a planet, either by a single en- 
counter or by r e peate dly perturbing the planet's orbit. 
IWeinberg et al.l (|1987t ) calculate that half-life of a long- 
period planet in the solar neighborhood due to the cu- 
mulative effect of many encounters with passing stars to 
be 

fAh\ ^ 10 4 AU ^ / 0T5MqPc^ \ 
S * 10Gyr {^) [—) { P ) ' (2) 

where a is the planet's semi-major axis. Assuming that 
the planet either formed in situ or migrated to its current 
location on a timescale much less than the system age, we 
can equate the timescale for a planet to become unbound 
to the age of the star, resulting in an estimate of the 
maximum distance at which a planet is likely to survive, 
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where is the age of the star, and provided that a ps <C 
at- Thus, the limiting semi-major axis for young disk 
stars will be set by the tidal limit, while the limiting semi- 
major axis for old field stars will be set by the cumulative 
effect of many stellar encounters. 

2.2. Stars in Clusters 

Mo st disk stars formed in small clusters (jThroop et al.l 
l2001f h Because the rate of stellar encounters depends 
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strongly on the local stellar density, the effects of stellar 
encounters are likely dominated by interactions while the 
star is young and still in a dense young stellar cluster. If 
planets form or migrate to long period orbits while the 
host star is still in a dense region of the cluster, then 
encounters with cluster members could place more strin- 
gent constraints than those discussed in H2.ll The im- 
portance of stellar encounters is linked to the dynamical 
evoluti on of the cluster in which the star formed. Fol- 
lowing lAdams et al.l (|2006[ ) , we separate the dynamical 
evolution of the cluster from the dynamics of close en- 
counters. 

2.2.1. Cross Sections for Impulsive Encounter with Single 

Stars 

First, we consider the ejection of planets in a sin- 
gle encounter. The timescale for a stellar encounter to 
eject a planet from its host star can be estimated by 
r -1 ~ (n+avoc), where is the local density of stars, 
Woo is the velocity at infinity of the passing star, a is 
the relevant cross section, and () indicates a time aver- 
age. For very long-period planets, encounters with single 
stars are impulsive when 
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where a is the semi-major axis, and Mi is the mass of the 
host star. In this regime, the cross section for ionizing 
the planet from the host star is 
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where M2 is the mass of the pas sing star , and G is the 
gravi tational constant (|Hutl il983:Eq. 13 of lFregeau et al.1 
12001 

In addition to planets being ejected, some planets will 
be exchanged into orbit around the passing star. For the 
sake of comparing to the cross section for ionization, we 
consider equal mass stars of mass M = Mi = M2, in 
which case 
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j!ut]li983;Eq. 14 of lFregeau etaD 12006) . For strongly 
impulsive encounters, the cross section for exchange is 
much less than the cross section for ionization, so ex- 
changes are relatively insignificant. 

2.2.2. Cross Section for Intermediate Encounter with 
Single Stars 

For planets at somewhat shorter orbital periods (a < 
1000 AU), the scattering behavior enters an alternative 
regime, where the cross sections for ionization and ex- 
change are comparable. Numerical simulations for a cir- 
cular planet encountering a single star find that 
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were M t = Mi + Mi + m is t he total mass (Eq. 15 and 
Fig. 3 of lFregeau et al.ll2006fh In this regime, when the 



planet is exchanged to the the passing star, the bind- 
ing energy of the planet before and after the star will be 
comparable (with a typical change of ~ 20%). There- 
fore, exchanges between stars of similar mass result in a 
planet orbit comparable to before the encounter. Among 
a coeval cluster of stars, such an exchange is unlikely to 
be recognizabl e, except in unusu al circumstances (e.g., 
PSR 1620+26: iFord et al.lf2000bh . 

2.2.3. Rate of Close Encounters for Stars in Clusters 

Simulations of star clusters can be used to calculate the 
rate of stellar flybys. Table 3 from lAdams et all (2006) 
provides empirical fits to the results of their simulations, 
with the form, 
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where b is the closest approach distance, and T an d 7 are 
fit parameters (see Table 3 of Adams et al. 2006). For a 
typical cluster of 100-1000 stars, they find that a typical 
star undergoes N(< 100 AU) ~ 0.2 — 2 or interactions 
with closest approach distance less than 100 AU over 
the course of a 10 Myr simulation. Similarly, they find 
N(< 1000 AU) — 1 - 14 encounters coming within 1000 
AU over the course of a 10 Myr simulation. 

2.2.4. Ionization Rate by Single Encounters with Single 

Stars 

We find the semi-major axis beyond which a planet will 
typically be ejected by a single encounter with a passing 
star by setting b = a-^j-K. Thus, a typical planet will 
be ejected if the encounter is impulsive (i.e., a > ai mp ) 
and the planet's orbit exceeds 
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where t is the time available for stellar encounters. Sim- 
ilarly, a typical planet will be ionized if t he encounters 
are in the intermediate regime identified bv lFregeau et all 
(2006) and 
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For planets that form concurrently with the star and 
in situ, the final factor of Eqs. (fT0|) & (fTTj) ranges from 
~ 0.5 — 6.9 for the 10 M yr simulations of s mall clus- 
ters with 100-1000 stars ([Adams et al.ll2006D . There- 
fore, planets that form (or migrate) beyond a,i mp during 
the earliest stages of the cluster would be ejected by a 
single impulsive stellar encounter. For some choices of 
cluster parameters (i.e., the "cold start" simulations of 
lAdams et al.l 120061 planets inward of ai mp but beyond 
a-ps,int could be ionized by passing stars before the clus- 
ter dissolves. 

The previous estimates assume that the planet forms in 
situ at the same time as the star. As we discussed in 
some theoretical models predict that planets might form 
close to a star and then be scattered outwards much later. 
For the same simulations, but considering planets that 
are placed at large separation 5 Myr into the cluster sim- 
ulations, the final factor of Eqs. (fT0|) & (jlip ranges from 
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~ 4.6 — 190 (co nsidering the subseq uent 5 Myr evolution 
of the clusters; lAdams et al.l (|2006f ). Therefore, planets 
that migrate to large separations af ter the cluster loses 
its gas (at 5 Myr in simulations of lAdams et al.|[2006f) 
could avoid being disrupted by stellar encounters at sep- 
arations of up to ~ max(ai mp , a P s,imp ) ~ 10 3 - 3 x 10 4 
AU. 

Most stars in a cluster are actually members of bi- 
nary (or higher order multiple) star systems. For the 
wide planetary orbits of interest, encounters with close 
binaries proceed similarly to encounters with single stars. 
The primary difference is that encounters which previ- 
ously resulted in the planet exchanging into orbit around 
the passing star now result in forming a triple system in- 
cluding the circumbinary planet. Encounters with binary 
stars that have larger orbital periods are more likely to 
result in ionization, as even exchange interactions often 
result in unstable triple systems that subsequently eject 
the lowest mass planet. Numerical simulations including 
a plausible distribution of binary mass ratios and peri- 
ods verify that the cross section for ionization by binary 
systems is roughly comparable to the sum of the cross 
sections for ionization by the individual stars. 

2.3. Summary of Survival of Long Period Planets 

In summary, interactions prior to the dissociation of 
the cluster can dominate the ionization of planets that 
form or migrate to large distances before the parent clus- 
ter dissociates. For planets that appear at large separa- 
tions near the birth of the cluster, survival for the life of 
the cluster can require a < 500 — 6000 AU, depending 
on the cluster properties. However, if planets are able to 
form or migrate to very large separations more than a 
few million years after the formation of the cluster, then 
they could survive at much larger distances. A delay of 
10 Myr is more than sufficient to allow planets to survive 
with semi-major axes approaching the tidal limit (or a ps 
for old star). 

In order to identify the most significant processes af- 
fecting planets in wide orbits, we compare several rele- 
vant timescales. The planet formation timescale is a key 
difference among the two leading models for planet for- 
mation: core ac cretion and gravitational instability (see 
lArmitagd [20071 for a recent summary). In the clas- 
sical core accretion model, accretion of gas onto a gi- 
ant planets at Jovian-like dista nces (~ 5 AU) lasts at 
least 8 Myr (|Pollack et al.lll996f) . However. lHaisch et all 
(|2001h estimates that the lifetimes of protostellar disks 
last at most w 3 Myr, and this timescale imposes a max- 
imum value on the core accretion timescale. This has 
inspired theoretical research to explore whether core ac- 
cretion may be able to act more rapidly. Indeed, varia- 
tions on the core accretion model are able to form plan - 
ets more rapidly. For example, iHubickvi et al.l (2005) 
find that Jupiter could have formed in 2.2 Myr, and 
iDodson-Robinson et al.l (|2008f ) present two new models 
where Saturn fo rms at ~ 12 AU in just under 3.5 Myr. 
iChambersI l[200fih envisages accretion times under 1 Myr. 
While further research is in order, it is unlikely that the 
core accretion timescale lies outside of 0.1 — 10 Myr. In 
the gravitational instability model, giant planets acquire 
their mass on much shorter timescales w hich are unlikely 
to surpass 10 4 yr (jDurisen et al.l [20071 and references 
therein) . 



After multiple planets have formed in a protoplane- 
tary disk, they can begin to interact with each other. 
Initially, the gas disk can suppress eccentricity excita- 
tion, but once the gas dissipates, the system can rapidly 
become orbit crossing. Alternatively, for systems con- 
taining planets that are more widely spaced, the planets 
can remain on low-eccentricity orbits for over a billion 
years before chaotic evolution leads to orbit crossing. In 
either case, once orbits cross, close encounters and eccen- 
tric orbits soon follow (|Chatteriee et al.ll2008[ ). Before a 
planet is to be scattered onto a very wide orbit or ejected 
from the system, it will undergo a series of many close 
encounters. The timescale depends on the planet masses 
and the orbital period of the more massive planet (that 
typically remains bound). For scattering by Jupiter, the 
median value of this timescale is ~ 10 4 yr, with a range 
that spans several magnitudes (see £)3.2j) . 

The process which limits the orbital separation of long- 
period planets will vary from system to system, depend- 
ing on the host star environment and the time until the 
planet is scattered onto a wide orbit. Most stars are be- 
lieved to form in star clusters, groups, or associations 
that dissociate on timescales of a few Myr. Simulations 
of star clusters show that most of the close encounters 
occur during the first few million years. If planets form 
in situ in wide orbits (e.g., by gravitational instability) 
or if they migrate or are scattered into wide orbits very 
rapidly, then planets beyond ~ 700 — 4000 AU could 
have been disrupte d by passing stars before the host 
cluster dissociates (jAdams etafl [20061 ). The width of 
the main sequence in open clusters suggests that stars 
in a larger cluster form within ~ 5 Myr of each other 
(jBonatto fc Bical l2009f ). Thus, it may be possible for 
planets formed via gravitational instability to survive at 
larger separations, if the planet orbits a star that formed 
a few Myr after the cluster's massive stars have started 
to remove the cluster's gas. For such planets and planets 
that are scattered into wide orbits shortly after the clus- 
ter dissociates, the tidal field of the molecular cloud can 
disrupt planets with separations exceeding ~ 10 4 — 10 5 
AU. For planets that are scattered into wide orbits well 
after the cluster has dissociated and the star has drifted 
away from the host molecular cloud, perturbations from 
passing stars and galactic tides can remo ve wide-period 
plane ts on timescales of ~ Gyr (Eq. [21 iHiguchi et al.l 
l2007h . 

For direct imaging to detect planetary systems like our 
solar system, it is important to focus on the closest young 
stars, where the inner working angle corresponds to a 
small physical separation. Most young stars in the imme- 
diate solar neighborhood are believed to have formed in 
small groups or associations, in which case the processes 
described in § 12.11 will not be relevant. However, these 
same direct imaging campaigns may overlook young gi- 
ant planets on wide orbits, as they lie beyond the detec- 
tor's outer working angle. When searching for planets at 
wide separations, it is possible or even advantageous to 
target more distant star forming regions. Thus, it would 
be possible to search for planets that formed in more 
substantial clusters that span a variety of ages. 

3. FORMATION OF LONG PERIOD PLANETS VIA 
PLANET-PLANET SCATTERING 
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When considering the survival of long-period planets 
in § 12. 1[ we applied arguments previously developed for 
the survival of Oort cloud comets. However, the mech- 
anisms responsible for creating the Oort cloud are un- 
likely to produce a large population of long-period plan- 
ets. Populating the Oort cloud is inefficient, as only 
~2% of the parent population scattered by Jupiter and 
Saturn end up in bound orbits with very long periods, 
rather than bei ng ejected from the system or co l liding 
with a planet (iDuncan et all 119871; iDones et all 12004 
iHiguchi et al.ll2006tlBrasser et alJl2007fT In our solar sys- 
tem, Uranus and Neptune may scatter planctcsimals into 
the Oort cloud more efficiently, largely due to their lower 
masses. However, they are not massive enough to scat- 
ter a Jupiter-mass planet into an Oort cloud-like orbit. 
Assuming that planet formation typically produces a few 
giant planets per star, a typical cluster might form a few 
systems where a giant planet is scattered outwards and 
be perturbed onto an Oort cloud-like orbit that will per- 
sist for billions of years. Although giant planets remain- 
ing in very long-period orbits for billions of years will 
be rare, the number of giant planets at large separation 
may be greatly increased around young stars, since giant 
planets that have been scattered outwards will typically 
spend less than 1 Myr on wide, eccentric orbits before 
becoming unbound. The chance of detecting such plan- 
ets is further enhanced since the direct detection of giant 
planets is most practical while they are still young and 
radiating the heat of accretion in the infrared. There- 
fore, we explore the possibility of detecting giant planets 
that have recently been scattered outwards, before they 
eventually become unbound. 

3.1. Methods 
3.1.1. Initial Conditions 

We investigate whether planet scattering may produce 
a population of detectable planets on wide orbits by per- 
forming n-body integrations of "dynamically active" (= 
will undergo instability) planetary systems that typically 
contain three giant planets and three less massive plan- 
ets. In the test particle limit of the less massive plan- 
ets, systems with three giant planets represent simple 
systems that can undergo a dynamical instability after 
an e xtended period of seemingly regu l ar orbital evolu- 
tion arzari fc Weidenschillind l2002t IChatteriee et al.l 
2008]). This is in contrast to systems containing only 
two planets, which have a sharp boundary dividing prov- 
ably stable configuration s from systems that rapidly re- 
sult in close encounters (|Gladmanl [i~993) . Three- planet 
systems have the additional advantage that they are 
more li kely to retain signatures of possible prior insta- 
bilities (|Veras fc Armitagel l2006f ). Although planetary 
systems may form even more planets, chaotic orbital 
evolution typically results in planets being ejected, ac- 
creting onto the star, or c olliding with each other until 
only 2 or 3 planets remain (|Papaloizou fc Terquem|[2001l : 
I Adams fc Laughlml 120031 ). Due to the strongly chaotic 
evolution of dynamically active multiple-planet systems, 
systems with th ree or more giant plane ts behave quali- 
tatively similar (jJuric fc Tre mainc 2008). 

We consider a realistic ensemble of post-formation ini- 
tial stellar masses, planet masses, and orbital separa- 
tions. We draw stellar masses from two different initial 



mass functions, and show that the results are broadly 
similar. We first utilize a simplistic initial mass function 
(IMF) with a power-law exponent of —1.35 with mass 
cut-offs at 0.3M Q and 3.0M Q , approximating the distri- 
bution of A through mid-M stars. We then apply rejec- 
tion s ampling in order to o btain a more realistic Miller- 
Scalo (|Miller fc Scaldfl979l) power-law distribution (with 
exponents of —0.25, —1.0 and —1.3 for the mass ranges 
O.3M - 1.0M Q , 1.0M© - 2.0M Q and 2.0M Q - 3.0M Q , 
respectively). We report all data and figures from the 
Miller- Scalo power-law distribution, and where compar- 
ison is appropriate, data for the simple power-law IMF 
in square brackets. The comparisons will show that the 
results are insensitive to this choice of IMF at the ~ 1% 
level. 

We randomly draw six planet masses from a distri- 
bution uniform in the log of the planet -star mass ratio 
as su ggested by radial velocity surveys (|Cumming et al.1 
2008). The upper limit (10 2 ) is based on the paucity of 
more massive planets around solar-type stars. The lower 
limit (9 x 10~ 6 ) is near a terrestrial mass and is well- 
beyond the reach of current direct imaging campaigns. 

Theoretically, the site of giant planet formation is be- 
lieved to occur predominantly beyond the location of 
the snow line. Therefore, we randomly sample the in- 
nermost planet's semi- major axis uniformly from 3 — 7 
AU. The lower bound roughly corresponds to the loca- 
tion of the snow line in the minimum mass solar neb- 
ula (|Lecar et al.l 120061 ). Various environments can yield 
snow li ne locations which vary over one order of mag - 
nitude (|Garaud fc L in 2007; [Kennedy fc Kenvoni r2008L 
Therefore, our range of semi-major axes for the inner- 
most planet allows for a variety of likely snow line loca- 
tions. 

We choose an initial spacing between neighboring plan- 
ets constant in units of the mutual hill radius, K . Our 
choice of K = 4.2 almost guranatees that each sys- 
tem will undergo dynamical instability, and does so 
wi thin a few million inn er orbital times (see Fig. 29 
of IChatteriee et aLll2008[ ). Thus, instabilities occur only 
after enough orbits that chaos has erased memory of the 
exact initial conditions. The value of K then effecively 
sets the timescale for instability in systems where we ex- 
pect close encounters, collisions and ejections to occur. 
Hence, our results should be interpreted in terms of the 
rate of detectable planets per dynamically active plane- 
tary system, rather than the rate relative to all stars or 
all planetary systems. Although our choice for the value 
of K sets the timescale for the initial onset of instability 
in our simulations, the final outcomes and the timescale 
for planets on wide, eccentric orbits to become unbound 
is not sensitive to K, but rather the mass of the inte- 
rior p lanet that eventually ejects the planet (|Pan fc Saril 
l20f)l . 

We randomly choose the orbital orientation and mean 
anomalies from uniform distributions, and assign the ini- 
tial eccentricities and inclinations also from uniform dis- 
tributions, bounded by < 0.01 and < 0.01°. Because we 
begin each system with three relatively closely spaced 
massive planets, the systems will typically self-excite 
larger eccentricities and inclinations. In this paper, we 
focus our attention on planets that have undergone re- 
peated close encounters with one or more planets before 
being scattered to very large (> 100 AU) separations. 
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Thus, our conclusions will not be sensitive to the exact 
choice of initial semi-major axes, eccentricities, and in- 
clinations. 

3.1.2. Integration Methods 

We integrate 1899 [2000] sets of initial conditions de- 
scribed above for 500 million years. We perform all in- 
tegrations with the Burlirsch-Stoer in tegrator from th e 
Mercury iV-body integration package (|Chambersl [1999) . 
In order to monitor the effects of ejected planets, plan- 
ets were removed from the system only upon reaching 
a semimajor axis of 10 6 AU from their parent star. All 
integrations included the effects of possible collisions of 
planets with one another and with the parent star, but 
did not include fragmentation. A different set of 2000 
integrations also included tidal effects from the galactic 
disk added to the fu nction mfo-us e r, acc ording to the 
prescription given bv lHiguchi et al.l (|2007f) : 



df 2 



-G M G) — 



(12) 



where vq — ^/AirGp is the vertical frequency, and p 
is the density in the solar neighborhood. We assumed 
p = O.lMgpc -3 (|Holmberg fc Flvnnl [2000h . In addition 
to the disk tide, there is a radial component of the galac- 
tic tide due to the galactic bulge and dependent on the 
Sun's angular speed around the galactic center. Since 
the radial galactic tide operates on a tim e scale an or- 
der of magnitude longer than the disk tide iHiguchi et al.1 
(|2007f ). we neglect this term. 

In Sj4l we will discuss the prospects for observing such 
planets. In this section, we focus on the dynamical prop- 
erties of these integrations. 

3.2. Results 

The numerical simulations described in [J3]result in sev- 
eral qualitative outcomes. In each case, at least one of 
the six planets is ejected, collides with its parent star, or 
collides with one another. Any of these outcomes causes 
the planet to be removed from the system. The follow- 
ing lists the percent of systems with a given number of 
bound planets remaining at the end of the simulations (at 
500 Myr), with the values in parentheses indicating the 
simulations which included the galactic tide: planets, 
1.4% [1.2% (1.2%)]; 1 planet, 15.4% [15.6% (15.5%)]; 2 
planets, 41.0% [41.3% (41.4%)]; 3 planets, 27.9% [27.9% 
(27.7%)]; 4 planets, 13.3% [13.0% (12.2%)]; 5 planets, 
1.16% [1.15% (1.95%)]; 6 planets, 0% [0% (0%)]. Hence, 
the broad dynamical effect of the galactic tide on these 
simulations, after 500 Myr and over distances of 10 6 AU, 
is negligible. Further, two-planet systems are the most 
likely outcome of our planet-planet scattering simula- 
tions, with three-planet systems being the second-most 
likely outcome; these results agree well with Fig. 6 of 
Adams & Laughlin (2003), who use 10 planets and a 
similar initial mass distribution. 

Although our simulations were run for 500 Myr, the 
vast majority of close encounters, collisions and ejections 
occur well before 100 Myr. The median time to the first 
close encounter between any two planets or planet is 
t ps 9.04 x 10 3 [9.15 x 10 3 ] yr, and the mean time to 
the first close encounter was t 1.33 x 10 6 [1.30 x 10 6 ] 
yr. We take one Hill Radii to be the maximum distance 
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Fig. 1. — Average number of planets vs. time for planets which: 
remain bound to their parent star (red/solid line), have been 
ejected (black/dotted), have collided with the star (blue/dashed 
line), and have collided with another planet (green/dot-dashed 
line). 



from a body which constitutes a close encounter. In or- 
der to help quantify the dynamical settling time of the 
systems, we plot the average number of planets remain- 
ing, planets ejected and planets suffering a collision as a 
function of time in Fig. [TJ The figure demonstrates that 
approximately half of the planets in each system still re- 
main after 2 x 10 7 yr, and that ejections (as opposed to 
collisions) become the predominant removal mechanism 
after 30 Myr. After 2 x 10 8 yr, on average, more plan- 
ets have been ejected than remain. Further, the fraction 
of planet-planet collisions are approximately equal to the 
fraction of planet-star collisions, but the latter dominates 
after 1 x 10 7 yr. At the end of our simulations, on av- 
erage we can expect a system to still have between 2-3 
planets. On even longer timescales we suspect that just 
a few more planets might be removed from the systems. 

Figure [3] presents snapshots of the semi-major axes, 
periastron distances, and apoastron distances for each of 
the bound planets remaining with masses greater than 
5 x 1O~ 4 M0 (a representative observational constraint) 
at times of 5, 10, 50, 100 and 500 Myr. Each row repre- 
sents a different system, and the five panels, in descend- 
ing order, contain 3583, 3476, 3227, 3164 and 3056 [3737, 
3623, 3364, 3298, and 3186] planets (from an initial pop- 
ulation of 12,000 planets). For the simple power law, 
without the mass constraint, the total number of bound 
planets at each of those 5 times would be [7537, 6805, 
5403, 5122 and 4783]. The horizontal lines indicate the 
periastron-apoastron range of each planet, and the points 
indicate the semi-major axes. Color indicates the initial 
ordering in terms of semi-major axis, as described in the 
caption. Because planets that have been scattered be- 
yond 100 AU are typically the least massive planets in 
the system, and there does not appear to be a significant 
correlation with the initial ordering (in terms of semi- 
major axes). By 5 and 10 Myr, the planet scattering 
has significantly excited planetary inclinations, but in- 
clinations of ss 90° < i < 150° are rarer than those of 
« 150° < i < 180°. We find 0.5% [0.5%] of all planets at 
5 Myr both harbor semi-major axes above 100 AU and 
are on retrograde orbits, but these planets are typically 
ejected after several 10 7 years. 

The planets of greatest interest to observational 
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Fig. 2. — Snapshots showing the orbital properties of the surviv- 
ing planets (of each system) with masses greater than 5 X 10 — 4 Mq 
at five times from 5 Myr (top) to 500 Myr (bottom). Within a 
panel, each row represents a different system. In the left pan- 
els, points indicate semimajor axes, and the horizontal lines in- 
dicate the pericenter and apocenter distance for each planet. In 
the right panels, points display inclination. The planet sym- 
bol colors indicate indicate initial proximity to the star, and are 
given by blue/circles (innermost), green/crosses, red/squares, ma- 
genta/triangles, salmon/diamonds and aqua/dots (outermost). 
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Fig. 3. — Cumulative histograms of the semimajor axis of 
the outermost bound planet at 5 Myr (black/solid lines), 10 
Myr (blue/dotted lines), 50 Myr (green/dashed lines), 100 Myr 
(red/dot-dashed lines) and 500 Myr (magenta/dot-dot-dot-dashed 
lines) . 



Fig. 4. — Cumulative histograms of the projected separation 
of the outermost bound planet at 5 Myr (black/solid lines), 10 
Myr (blue/dotted lines), 50 Myr (green/dashed lines), 100 Myr 
(rcd/dot-dashed lines) and 500 Myr (magenta/dot-dot-dot-dashed 
lines) . 

searches are those on "wide" orbits, with semimajor axes 
of tens, hundreds or thousands of AU. The number of 
planets surviving beyond 1000 AU declines rapidly from 
10 Myr to 50 Myr, but even at 500 Myr, planets re- 
main beyond 1000 AU. In Fig. [3l we present a cumu- 
lative distribution for the semimajor axis of the outer- 
most surviving planet from 10 1 ' 5 — 10 7 AU at 5 Myr 
(black/solid lines), 10 Myr (blue/dotted lines), 50 Myr 
(green/dashed lines), 100 Myr (red/dot-dashed lines) 
and 500 Myr (magenta/dot-dot-dot-dashed lines). The 
ejection of many planets (especially between 10 Myr and 
50 Myr) partially accounts for the difference between the 
two sets of curves (compare the 5 and 10 Myr curves to 
the 50, 100 and 500 Myr curves). Fewer planets reside 
beyond 40 AU at 5 Myr than at 10 Myr, because many 
systems are still undergoing significant dynamical evo- 
lution at 10 Myr. Between 10 Myr and 100 Myr, the 
number of planets on wide orbits decreases, since the 
rate of planets being scattered outwards into wide orbits 
is less than the rate of planets already on wide orbits be- 
coming unbound. At 100 Myr, less than 0.1% of systems 
harbor planets with semi-major axes beyond 10 5 AU. Al- 
though small, such planets could be present in a cluster 
of thousands of stars. 

Unfortunately, determining the semi-major axes of 
planets on very wide orbits is difficult because the orbital 
timescale greatly exceeds the typical time span of obser- 
vations. Instead, direct imaging measures the star-planet 
separation projected onto the plane of the sky. In order 
to evaluate the detectability of such planets, we present 
cumulative distributions for the projected separation in 
Fig. 3J In order to improve our statistics, we calcu- 
late the projected separation for each of five times and 
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Fig. 5. — Histogram of the eccentricity of the outermost bound 
planet at 5 Myr (black/solid lines), 10 Myr (blue/dotted lines), 
50 Myr (green/dashed lines), 100 Myr (red/dot-dashed lines) and 
500 Myr (magenta/dot-dot-dot-dashed lines). Bound planets with 
nearly parabolic orbits are prevalent for several tens of Myr, but 
eventually become unbound. 

ten viewing geometries (by drawing the mean anomaly, 
pericenter angle, and cos(z) from uniform random distri- 
butions). For nearly circular orbits, the planet-star sep- 
aration remains close to the semi-major axis throughout 
the orbit, and the projected separation averages to be 
~ 0.7 times the semi-major axis. For highly eccentric 
orbits, the planet spends most of its time near apocen- 
ter, so that the average projected separation can exceed 
the semi-major axis, depending on the orientation of the 
orbit relative to the line of sight. We find that the two 
effects approximately cancel for planets on wide orbits. 
The cumulative distributions in Figs. [3] and 0] differ by 
just a few percent. 

Next, we investigate the eccentricities of the planets 
that have been scattered onto wide orbits. Previous 
studies have compared the eccentricities of planets af- 
ter planet scattering to radial velocity observations. For 
example, two-planet scattering rarely results in an inner 
planet eccentricity grea ter than 0.8 followi ng the ejec- 
tion of another planet (|Ford &: Rasiol l2~008h . Similarly, 
three-planet scattering results in only a modest frac- 
tion of inner planets w ith eccentricities greater than 0.8 
([Chatter jee et al.ll2008r ). In this study, we focus instead 
on the eccentricities of planets that have been scattered 
outwards. Although many of these planets will eventu- 
ally become unbound, we consider the eccentricity dis- 
tribution of all planets still bound at each of several ages 
in Fig. The black/solid, blue/dotted, green/dashed, 
rcd/dot-dashcd and magenta/dot-dot-dot-dashed lines 
on this histogram correspond to 5, 10, 50, 100 and 500 
Myr. The last bin demonstrates that the number of plan- 
ets on nearly parabolic (e > 0.95) orbits decreases with 



time, where the greatest difference is between the 10 and 
50 Myr curves. Taken together, Figs. [3] and [5] indicate 
that planets on the verge of ejection typically have nearly 
linear orbits, spending the vast majority of the time near 
apocenter, but receiving a "kick" once each orbit when 
they bri efly sweep into the inner portion of the planetary 
system dPanfc Saril feOOT). In our simulations, less mas- 
sive planets are preferentially scattered into wide orbits. 
The majority of these planets are less than a Jupiter 
mass, but only a few percent exceed AMj up . 

Figures [3][5] do not include planets that have already 
been ejected onto hyperbolic orbits (e > 1). An un- 
bound planet will quickly become undetectable to a coro- 
nagraphic direct imaging search with a relatively small 
field of view. Assuming a planet with an ejection veloc- 
ity of ~ GM+/a, the ejected planet would have traveled 
more than 10 AU away from the host star (beyond at) 
in ~ 10 5 years. As a result, directly detecting an ejected 
planet still in the vicinity of its parent star is improb- 
able. However, such planets may still be detectable as 
free-floating planets by direct imaging of more distant 
star forming regions which image large fields of view. 

4. DETECTABILITY OF LONG PERIOD PLANETS 

Here, we investigate the plausibility of detecting plan- 
ets at large separations with direct detection campaigns. 
We take snapshots of each planetary system from the 
previous section at sev eral times: 5, 10, 50, 100 and 500 
Myr, all times for which Baraffc et alJ (|2008h provide the- 
oretical models. The first two snapshots are particularly 
relevant for direct imaging of very young stars, typically 
of star forming regions ~ 100 pc from the Sun. In those 
locations, models predict giants planets will be signifi- 
cantly more luminous in J band than a blackbody of the 
same temperature. The 100 and 500 million year snap- 
shots are particularly relevant for observations of nearby 
stars. At these later ages, models predict that planets 
have cooled enough that the J band luminosity is dra- 
matically reduced, so L band observations may become 
more favorable. 

For direct imaging to detect planetary systems like our 
solar system, it is important to focus on the closest young 
stars (at w 100 pc), where the inner working angle corre- 
sponds to a small physical separation. Most young stars 
in the immediate solar neighborhood are believed to have 
formed in small groups or associations, in which case the 
processes described in $2.2\ will not be relevant. How- 
ever, these same direct imaging campaigns may overlook 
young giant planets on wide orbits, as they lie beyond the 
detector's outer working angle. When searching for plan- 
ets at wide separations, it is possible or even advanta- 
geous to target more distant star forming regions. Thus, 
it would be possible to search for planets that formed in 
more substantial clusters at a variety of ages. 

4.1. Detection Criteria 

Snapshots of each planetary system from the previous 
section at 10, 50, 100 and 500 Myr, are plotted in terms 
of contrast and angular separation in Fig. [6l Contrast is 
calculated by comparing the absolute magnitude of the 
plane t to that of its host. We use the atmospheric mod- 
els of iGirardi etldl (|2002T ) and lBaraffe etal] (pOOl ) (and 
Baraffe et al. 2008, private communication) to find the 
brightness of the stars and exoplanets respectively. In 
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Fig. 6. — Contrast versus angular separation for a population of exoplanets as a function of age. Data points are color-coded/grayscaled 
according to eccentricity. Circles indicate that all six planets in a given system remained gravitationally bound. Triangles indicate that at 
least one of the planets became unbound. Dashed lines indicate the approximate range of detectability based on current and ground-based 
direct imaging searches. Companions located exterior to the FOV edge are not detectable. The estimates for GPI and the TMT PFI arc 
somewhat conservative, as they could achieve deeper effective contrasts by a factor of order unity. 



Fig. [51 we assume that a cluster of stars, each with met- 
alicity Z = 0.019, is located at an average distance of 
100 pes from the Sun with an RMS scatter of 15 pes. 
In Fig. we place the same stars at 60 pes with an 
RMS scatter of 20 pes. In each case, their positions 
within these limits are randomized. For reference, we 
also include approximate detection se nsitivities from re- 
cent AO-coronagraph surveys (namely iBiller et al.ll2007l; 
Kasper et al.l 120(371 : lLafreniere et al.l 120071 : iNielsen et all 
2007F as well as future instruments, such as G PI - The 
Gemini Planet Imager (Macinto sh et al.ll2006af) and the 
proposed TMT PFI - Thirty Meter Teles cope Planet For- 
mation Imager (|Macintosh et al.ll2006bl ). 

Models predict that giant planets will be several orders 
of magnitude more luminous in the near-IR bands com- 



pared to a black-body of the same temperature at such 
young ages. The first two snapshots (10+50 Myrs, top 2 
panels of Fig. [6]) are particularly relevant for searching 
members of star forming regions. At these two times we 
selected the J-band (A = 1.25 /im) because it offers the 
brightest planets, best spatial resolution, and faintest sky 
background in a regime where the atmosphere can still 
be adequately controlled with AO. 

The 100 and 500 Myrs snapshots are more relevant for 
observations of nearby systems, such as Hyades or the 
Ursa Majoris cluster. At these late stages, models predict 
that planets have cooled enough that the J-band luminos- 
ity is dramatically reduced. In this case, L-band observa- 
tions, A = 3.8 fJ XOj may beco me favorable for large angu- 
lar separations (Hcinze 2007). There is, however, a trade- 
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Fig. 7. — Contrast versus angular separation for a nearby population of exoplanets in the J- (upper panels) and L (lower panels)-b&nds. 
Planets are brighter in the L-band at older ages, but the sky background is at least 10 — 6 times as bright as the brightest stars after 
subtraction and places a hard-limit on ground-based sensitivities. GPI is unlikely to conduct observations in the L-band. 



off, because the L-band has the disadvant age of a bright 
and rapidly fluctuating sky background l|Phillips et alj 
1999), which ultimately limits sensitivity (Fig. [7]). Chop- 
ping and nodding techniques can only subtract out this 
source of noise (plus instrument and detector noise) to 
one part in rs 10 5 , even when operating at the photon 
noise limit. In the best case scenario, this translates to 
contrast levels of w 10 -6 , i.e. when targeting stars with 
apparent magnitudes L < 4. The prospects are better 
by nearly an order of magnitude when conducting obser- 
vations at the South Pole, but constructing single dish 
telescopes with sufficient angular resolution is currently 
difficult given the harsh conditions. 

It is also worth noting that the sensitivity of special- 
ized high-contrast imaging instruments is non-monotonic 
with angular separation. This is due to the explicit re- 



moval of instrument scattered light within the control- 
lable spatial frequencies of the AO system deformable 
mirror. The number of actuators and their spacing de- 
fines this primary search area - the so-called "dark-hole" 
(|Trauger fc Traubl [2001 . GPI's field of view (FOV) is 
3.6". Current high res olution instruments have some- 
what larger FOV's (e.g. lHavward et al.|[200ll) . 

Finally, in Fig. [SJ we plot the fraction of orbit realiza- 
tions of all planets which would be detectable according 
to Figs. [6] and [7] as a function of age, band, and sepa- 
ration. We emphasize that this fraction is representative 
of all planets in dynamically active systems (those which 
undergo instability), according to our initial conditions. 
Current technology should allow for the direct imaging 
detection of planets at ages < 100 Myr, whereas GPI 
and TMT PFI are predicted to discover several times 
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Fig. 8. — Percent of all planets from an enhanced distribution 
(where 10 viewing geometries for each initially 6-planet system are 
included) which are detectable within the bounds imposed in Figs. 
[6] and [7] according to current technology (solid/black curve), GPI 
(blue/dotted curve) and TMT PFI (green/dashed curve) in the J- 
band for a cluster located at an average distance of 100 pes from 
the Sun with an RMS scatter of 15 pes. All isolated symbols re- 
fer to a cluster located at an average distance of 60 pes from the 
Sun with an RMS scatter of 20 pes and represent: J-band, cur- 
rently detectable (red crosses); J-band, GPI detectable (magenta 
squares); J-band TMT PFI detectable (salmon triangles); L-band 
currently detectable (aqua diamonds); and L-band TMT PFI de- 
tectable (grey asterisks). 

more planets. Both instruments are likely to discover 
the greatest number of planets at an age of ss 50 Myr. 
Importantly, current J-band (red crosses) and L-band 
(aqua diamonds) observations demonstrate the slight im- 
provement the latter affords at later times; the projected 
TMT PFI observations in the J-band (salmon triangles) 
and L-band (grey asterisks) demonstrate the marked im- 
provement the L-band provides at both 100 Myr and 500 
Myr. 

5. DISCUSSION 

We have shown that planets could survive in very- wide 
orbits. For the relatively young stars typically surveyed 
by direct imaging campaigns, planets could survive at 
distances of up to the tidal limit. As for a potential 
formation mechanism, we consider the role of planet- 
planet scattering, with a focus on the implications for 
long-period planets. Because the process of planet-planet 
scattering typically extends for tens of millions of years, 
most stars will not remain in a high-density stellar envi- 
ronment long enough for close stellar encounters to be ef- 
fective at removing long-period planets (with a < at, a ps ; 
see Eqs. [Hand [3]). 

We performed direct n-body simulations of dynami- 
cally active planetary systems, including the effects of 
galactic tides, and show the planet-scattering can in- 
deed form long-period planets. Of particular interest, 



we found that at ages of ~ 10 — 50 Myr, there is a siz- 
able population of planets on wide and highly eccentric 
orbits. Although most of these planets will eventually be 
ejected, they can persist on wide orbits for several tens of 
millions of years. We caution that our models predict a 
significant temporal evolution of this population. There- 
fore, statistical analyses of direct imaging surveys should 
account for the possibility of a time-variable population 
of planets on wide orbits. 

We considered the detectability of the planets scattered 
into long-period orbits in our simulations. Planet models 
are highly uncertain, largely due to the uncertainty in 
initial conditions (e.g. lMartin et al.l|200 7i). Nevertheless, 
we show that even assuming "cold-start" models, direct 
imaging surveys are already sensitive to giant planets 
in nearby star forming regions. Future direct imaging 
surveys will become sensitive to both closer planets and 
planets with smaller luminosities. 

Recen t interest in the can didate three-planet HR 8799 
system (|Marois et a l. 2008) motivates us to consider the 
multiplicity of planets from the same system all in bound 
wide orbits. The properties (such as eccentricity and in- 
clination) of the orbits of those three planets are largely 
unconstrained at the present time. The planet masses 
are approximated to lie in the 5 - l3Mj up range and the 
projected separations are estimated to be 24, 38, and 68 
AU (Marois et al. 2008). In our simulations of §3, the 
initial semi-major axis of the outermost planet is typi- 
cally 15-30 AU, so any planets beyond 40 AU have un- 
dergone significant scattering. In order to evaluate the 
frequency of multiple planets at wide separations due to 
planet scattering, we tally how many systems contain 
multiple planets on bound orbits with semi-major axes 
exceeding 40 AU and masses exceeding 5 x 1O~ 4 M . At 
(5, 10, 50, 100, 500) Myr, (71, 53, 30, 24, 18) [76, 59, 
34, 26, 19] systems contain two such planets and (5, 1, 
1, 1, 0) [7, 1, 1, 1, 0] systems contain three or more such 
planets (see FiglHJ). 

These results suggest that 3.8% (0.9%) of dynamically 
active planetary systems could contain multiple massive 
planets on wide orbits at ages of 5 Myr (500 Myr and be- 
yond). Further, 66% (40%) of systems containing multi- 
ple planets on wide orbits with ages less than 10 Myr (100 
Myr) will become unstable on longer timescales. This 
model predicts that multiple planets on wide orbits will 
become less common with increasing stellar age (Fig. [9]). 
Figure fTOl also demonstrates this trend by displaying the 
percent of systems with each adjacent pair of the out- 
ermost 3 planets having a projected separation within 
factors appearing on the X-axis. Black, blue, green, red 
and magenta curves refer to data at 30, 40, 50, 100, and 
150 Myr, ages chosen to corres pond to the likely age of 
the three planets in HR 8799 (|Marois et aTl|2008l ). On 
this plot, < 1% of all our systems simulated would be as 
"packed" as HR 8799 (with an X value of 38/24 = 1.58), 
suggesting that the planets in that system are unlikely 
to have been scattered to those locations wi thout reso- 
nant locking (|Fabryckv fc Murrav-Cla v 2008) or migrat- 
ing through a disk. 

Young systems containing multiple planets at wide sep- 
arations would be expected to have sizable eccentricities, 
even larger than average of extrasolar planets. With 
time, 2-planet systems with large eccentricities tend to 
be destroyed, so the average eccentricity decreases to a 
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Fig. 9. — Statistics of systems with multiple bound planets 
all on wide orbits (projected separation > 40 AU) as a function 
of age. The upper panel plots the fraction of all systems with two 
(black/solid line and dots), three (blue/dotted line and crosses) and 
four (green/dashed line and squares) planets remaining beyond 40 
AU. The red/dot-dashed line with triangles and the magenta/dot- 
dot-dot-dashed line with diamonds refer to systems with two or 
three massive (> 0.3Afj) planets remaining. The lower panel show- 
cases the median eccentricities (dots and black/solid lines) and sine 
of the relative inclinations (crosses and blue/dashed lines) of all 
systems with two planets on wide orbits, and their 25-75 percentile 
ranges. The blue/dashed lines are slightly offset from the other 
lines for clarity. 



more modest value (0.14), very similar to that observed 
for multi ple planet systems discovered by radial velocity 
searches (Wri ght et al.| [2008) . The average relative incli- 
nation between two planets on wide orbits also decreases 
with time from 28.4% to 9.5% degrees between 10 and 
500 Myr. These results can be used to inform the inter- 
pretation and dynamical mode ling of multiple planet sys- 
tems found by direct imaging (|Fabryckv fc Murrav-Clavl 
12008ft . 

In the future, comparing the frequency of long-period 
planets as a function of stellar age could provide a sig- 
nificant constraint on planet formation models. Such ob- 
servations could determine the timescale for the onset of 
strong instabilities in multiple planet systems and place 
constraints on the typical number of giant planets formed 
(which may be greater than the number of giant planets 
around mature stars). Such constraints would be par- 
ticularly interesting, when combined with constraints on 
free-floating planets from direct imaging of star forming 
regions and/or gravitational microlcnsing observations. 

In our simulations excluding the galactic tide, 40% 
[42%] of all planets (totaling 13.6% [13.7%] of the plan- 
etary mass) were ejected. Of the planets ejected, 21.2% 
(28.9%) [21.0% (28.9%)] had a mass (mass ratio) greater 
than O.MIjup (0.3M Jup /M G ) and 54.1% (60.1%) [54.2% 
(60.1%)] of planets had a mass (mass ratio) between 5M ffi 
(5M e /M Q ) and 0.3M Jup (0.3M Jup /M Q ). The simula- 
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Separation Factor 

Fig. 10. — Mutual separations of outermost 3 planets in each 
system, such that each adjacent pair of the outermost 3 planets 
has a projected separation within a factor of the values on the X- 
axis. Black/solid, blue/dotted, green/dashed, red/dot-dashed and 
magenta/dot-dot-dot-dashed lines refer to snapshots at 30, 40, 50, 
100, and 150 Myr. 



tions which included the effects from the galactic tide 
yielded nearly identical results, varying the above values 
by at most by a few percent. All these ejected plan- 
ets would become free-floating planets. After they fade 
and become undetectable by direct imaging, the best 
chance for detecting this population of planets is via 
gravitational microlcnsing. Although the optical depth 
for microlcnsing is independent of the mass function of 
free-floating planets, the event rate and event timescale 
are sensitive to the mass function. The event rate for 
ejected microlensing planets which are more massive 
than 0.3M Jup (0.3M Jup /M Q ) is 56.6% (64.4%) [56.2% 
(64.4%)] of the total event rate for all unbound masses, 
but other mass bins will produce different event rates. In 
order to better quantify this dependence on mass, we bin 
the masses uniformly in log M and for sufficiently high 
sampling, estimate that a power-law index of [« 0.37] 
approximates the positive correlation between event rate 
(as the dependent variable) and mass. 

In reality, there would be additional microlensing 
events due to 1) planets still bound to their parent stars, 
and 2) small planets or protoplanets there were not in- 
cluded in our simulations. Therefore, the above numbers 
should be interpreted as a lower limit. Regarding the 
contribution from bound planets, 38.0% [39.9%] of all the 
initial planets (totaling 81.5% [81.3%] of the planetary 
mass) remain bound. Of these, 76.5% (84.1%) [76.2% 
(84.1%)] had a mass (mass ratio) greater than 0.3Mj up 
(O.UIjup/Mq), and accounted for 93.8% (96.4%) [93.6% 
(96.4%)] of the overall microlcnsing event rate. Unfor- 
tunately, existing surveys are prone to missing short- 
duration events and hence biased towards finding more 
massive planets (especially for free-floating planets where 
there is not a longer duration stellar microlensing event 
to trigger intensive observations). These results suggest 
that a microlensing survey that efficiently detected the 
short events due to free-floating terrestrial planets would 
provide a useful diagnostic for the rate at which terres- 
trial planets form in dynamically active systems (relative 
to the rate of giant planet formation in similar systems). 
When combined with statistical information about the 
frequency of terrestrial and giant planets in mature plan- 
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etary systems, this synergy would provide a probe of the 
chaotic phase of planet formation and a test of the planet- 
planet scattering model for exciting eccentricities of ex- 
trasolar planets. 

We considered the detectability of the planets scat- 
tered into long-period orbits in our simulations. Indeed, 
planet models are highly uncertain largely due to th e 
uncertainty in initial conditions fe.g. lMartin et al1l2007t ). 
We show that even assuming "cold-start" models, direct 
imaging surveys are already sensitive to giant planets 
in nearby star forming regions. Future direct imaging 
surveys will become sensitive to both closer planets and 
planets with smaller luminosities. 
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Sigurdsson, Jonathan Tan, and Scott Tremaine for shar- 
ing their insights. The authors acknowledge the Uni- 
versity of Florida High-Performance Computing Center 
for providing computational resources and support that 
have contributed to the research results reported within 
this paper. 



REFERENCES 



Adams, F. C, Proszkow, E. M., Fatuzzo, M., & Myers, P. C. 2006, 

ApJ, 641, 504 
Adams, F. C, & Laughlin. G. 2003. Icarus. 163, 290 
Armitage, P. J. 2007, arXiv:astro-ph/0701485 

BarafTc, I., Chabrier, G., Barman, 1. s., Aiiard, F., & Hauschildt, 

P. H. 2003, A&A, 402, 701 
BarafTc, I., Chabrier, G., & Barman, T. 2008, A&A, 482, 315 
Biller, B. A., et al. 2007, ApJS, 173. 143 
Bonatto, C, & Bica, E. 2009, arXiv:0901.0833 
Boss, A. P. 2006, ApJL, 637, L137 

Brasser, R., Duncan, M. J., & Levison, H. F. 2007, Icarus, 191, 413 
Burrows, A., et al. 1997, ApJ, 491, 856 

Cai, K., Durisen, R. H., Michael, S., Boley, A. C, Meji'a, A. C, 
Pickett, M. K., & D'Alessio, P. 2005, Protostars and Planets V, 
8155 

Chambers, J. E. 1999, MNRAS, 304, 793 
Chambers, J. E. 2006, ApJ, 652, L133 

Chatterjee, S., Ford, E. B., Matsumura, S., & Rasio, F. A. 2008, 
ApJ, 686, 580. 

Cumming, A., Butler, R. P., Marcy, G. W., Vogt, S. S., Wright, 

J. T., & Fischer, D. A. 2008, PASP, 120, 531 
de la Fuente Marcos, R., & de la Fuente Marcos, C. 2008, ApJ, 

672 342 

Debes', J. H., & Sigurdsson, S. 2006, A&A, 451, 351 
Dodson-Robinson, S. E., Bodenheimer, P., Laughlin, G., Willacy, 

K., Turner, N. J., & Beichman, C. A. 2008, ApJ, 688, L99 
Dones, L., Weissman, P. R., Levison, H. F., & Duncan, M. J. 2004, 

Comets II, 153 

Duncan, M., Quinn, T., & Tremaine, S. 1987, AJ, 94, 1330 
Durisen, R. H., Boss, A. P., Mayer, L., Nelson, A. F., Quinn, T., 

& Rice, W. K. M. 2007, Protostars and Planets V, 607 
Fabrycky, D. C, & Murray-Clay, R. A. 2008, arXiv:0812.0011 
Ford, E. B., Havlickova, M., & Rasio, F. A. 2001, Icarus, 150, 303 
Ford, E. B., Joshi, K. J., Rasio, F. A., & Zbarsky, B. 2000a, ApJ, 

528, 336 

Ford, E. B., Kozinsky, B., & Rasio, F. A. 2000b, ApJ, 535, 385 
Ford, E.B., & Chiang, E.I. 2007, ApJ, 661, 602. 
Ford, E. B., Rasio, F. A., & Yu, K. 2008, ApJ, 686, 621. 
Fregeau, J. M., Chatterjee, S., & Rasio, F. A. 2006, ApJ, 640, 1086 
Fregeau, J. M., Cheung, P., Portegies Zwart, S. F., & Rasio, F. A. 

2004, MNRAS, 352, 1 
Garaud, P., & Lin, D. N. C. 2007, ApJ, 654, 606 
Girardi, L., Bcrtelli, G., Bressan, A., Chiosi, C, Groenewegen, 

M. A. T., Marigo, P., Salasnich, B., & Weiss, A. 2002, A&A, 

391, 195 

Gladman, B. 1993, Icarus, 106, 247 
Goldreich, P., et al. 2004, ARA&A, 42, 529. 

Gomes, R., Levison, H. F., Tsiganis, K., & Morbidelli, A. 2005, 
Nature, 435, 466 

Haisch, K. E., Jr., Lada, E. A., & Lada, C. J. 2001, ApJ, 553, L153 
Hayward, T. L., Brandl, B., Pirger, B., Blacken, C, Gull, G. E., 

Schoenwald, J., & Houck, J. R. 2001, PASP, 113, 105 
Heinze, A. N. 2007, Ph.D. Thesis, 

Higuchi, A., Kokubo, E., & Mukai, T. 2006, AJ, 131, 1119 
Higuchi, A., Kokubo, E., Kinoshita, FL, & Mukai, T. 2007, AJ, 134, 
1693 

Holmberg, J., & Flynn, C. 2000, MNRAS, 313, 209 
Hubickyj, O., Bodenheimer, P., & Lissauer, J. J. 2005, Icarus, 179, 
415 



Hut, P. 1983, ApJ, 268, 342 

Juric, M, & Tremaine, S. 2008, ApJ, 686, 603. 

Kalas, P. et al. 2008 Science Express DOI: 10. 1126/science. 1166609 
Kasper, M., Apai, D., Janson, M., & Brandner, W. 2007, A&A, 
472, 321 

Kennedy, G. M., & Kenyon, S. J. 2008, ApJ, 673, 502 
Lafrcniere, D., et al. 2007, ApJ, 670, 1367 
Laughlin, G., & Adams, F. C. 1998, ApJ, 508, L171 

Lecar, M., Podolak, M., Sasselov, D., & Chiang, E. 2006, ApJ, 640, 
1115 

Levison, H. F., Lissauer, J. J., & Duncan, M. J. 1998, AJ, 116, 
1998 

Macintosh, B., et al. 2006a, SPIE, 6272, 62720L 
Macintosh, B., et al. 2006b, SPIE, 6272, 62720N 
Mar ley, M. S., Fortney, J. J., Hubickyj, O., Bodenheimer, P., & 

Lissauer, J. J. 2007, ApJ, 655, 541 
Marois, C. et al. 2008 Science Express 

DOI:10.1126/science.ll66585 
Martin, R. G., Lubow, S. H., Pringle, J. E., & Wyatt, M. C. 2007, 

MNRAS, 378, 1589 
Marzari, F., & Weidenschilling, S. J. 2002, Icarus, 156, 570 
McCarthy, C, & Zuckerman, B. 2004, AJ, 127, 2871 
Meji'a, A. C, Durisen, R. H., Pickett, M. K., & Cai, K. 2005, ApJ, 

619, 1098 

Miller, G. E., & Scalo, J. M. 1979, ApJS, 41, 513 
Moorhead, A. V., & Adams, F. C. 2005, Icarus, 178, 517 
Morbidelli, A., Levison, H. F., Tsiganis, K., & Gomes, R. 2005, 
Nature, 435, 462 

Nielsen, E. L., Close, L. M., Bille r. B. A.. M asciadri. E., & Lenzen, 

R. 2007, ArXiv e-prints, 706, arXiv:0706.4331 
Pan, M., & Sari, R. 2004, AJ, 128, 1418 

Papaloizou, J. C. B., & Terquem, C. 2001, MNRAS, 325, 221 
Phillips, A., Burton, M. G., Ashley, M. C. B., Storey, J. W. V., 

Lloyd, J. P., Harper, D. A., & Bally, J. 1999, ApJ, 527, 1009 
Pollack, J. B., Hubickyj, O., Bodenheimer, P., Lissauer, J. J., 

Podolak, M., & Greenzweig, Y. 1996, Icarus, 124, 62 
Rafikov, R. R. 2005, ApJL, 621, L69 
Rafikov, R. R. 2007, ApJ, 662, 642 

Rasio, F. A., & Ford, E. B. 1996, Science, 274, 954 

Scharf, C. & Menou, K. 2008, submitted to ApJL. arXiv:0811.1981 

Stamatellos, D., & Whitworth, A. P. 2008, arXiv:0809.5042 

Thomess, E.W., Duncan, M.J., Levison, H.F. 1999, Nature402, 635. 

Thomess, E.W., Duncan, M.J., Levison, H.F. 2002, AJ123, 2862. 

Throop, H.B. et al. 2008, Science 292, 1686. 

Trauger, J. T., & Traub, W. A. 2007, Nature, 446, 771 

Tremaine, S. 1993, Planets Around Pulsars, 36, 335 

Tsiganis, K., GOmes, R., Morbidelli, A., Levison, H.F. 2005 

Nature435, 459. 
Veras, D., & Armitage, P. J. 2004, MNRAS, 347, 613 
Veras, D., & Armitage, P. J. 2006, ApJ, 645, 1509 
Weidenschilling, S. J., & Marzari, F. 1996, Nature, 384, 619 
Weinberg, M. D., Shapiro, S. L., & Wasserman, I. 1987, ApJ, 312, 

367 

Wright, J. T., Upadhyay, S., M arcy, G. W., Fisc her. D. A., Ford, 
E. B., & Johnson, J. A. 2008, arXiv:0812.1582 



